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CN ' Abstract 

The small x behavior of the flavor non-singlet gi structure function is 
' analysed numerically by taking into account the all-order resummation of 

as In^ X terms. We include a part of the next-to-leading logarithmic cor- 



rections coming from the resummed "coefficient function" which are not 
considered in the literatures to respect the factorization scheme indepen- 
dence. The resummed coefficient function turns out to give unexpectedly 
. large suppression effects over the experimentally accessible range of x and 

Q^. This fact implies that the higher order logarithmic corrections are very 
important for gi in the small x region. 



X 



1 Introduction 

The behavior of the structure function in the small Bjorken x region has recently 
received much attention of the physicists |l|]. This region corresponds to the so- 
called Regge limit. So we naively expect that the soft physics (Regge theory) 
may explain the small x behavior of the structure function. However the steep 
rise of the structure function in this region observed by the HERA experiments 
contradicts with this naive expectation. The physics at small x is now one of 
the most interesting subjects and many people believe that this problem could 
be handled in the context of the QCD perturbation theory 0. 

In the case of the polarized structure function gi, we have not yet had data 
at very small x . However it is very important and desirable to know the small 

*Talk presented by J. Kodaira. To appear in: Proceedings of the Workshop on Deep Inelastic 
Scattering off Polarized Targets: Theory Meets Experiment, DESY Zeuthen, September 1-5, 
1997. 
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X behavior of gi in the hght of the Bjorken and EUis-Jaffe sum rules. The naive 
Regge theory tells us, 

^1 ~ x", < a < 0.5 , 

at small x and and this implies a "flat" input parton density. Although 
the recent detailed analyses @] of the existing data in the DGLAP approach 
do not necessarily require an input function which contradicts with the naive 
Regge prediction, it is inevitable to take into account the effects from the Inx 
corrections, which appear in the perturbative calculations, to get more precise 
predictions on the small x behavior of gi [^]. 

Bartels, Ermolaev and Ryskin ^ have given the resummed expression for the 
gi structure function by using the Infra-Red Evolution Equation and confirmed 
the old result by Kirschner and Lipatov 0. They claim that the resummation 
effects may be very important. On the other hand, the numerical analysis by 
Bliimlein and Vogt shows that there are no significant contributions from the 
resummation of the leading logarithmic (LL) corrections in the HERA kinemat- 
ical region ( x ~ 10~^). This controversial aspect might come from the fact that 
the resummed "coefficient function" is considered in Ref. but not in Ref. f^. 
Bliimlein and Vogt did not include the resummed coefficient function because 
it falls in the next-to-leading logarithmic (NLL) corrections and depends on the 
factorization scheme adopted. It is also to be noted that the evolution, in gen- 
eral, strongly depends on the input parton densities. If one chooses a steep input 
function, the perturbative contribution will be completely washed away. So it 
will be interesting to see also the sensitivity of the results to the choice of the 
input densities. 

In this report, we discuss the numerical impact of the all-order Inx resumma- 
tion on the small x behavior of gi (non-singlet part). We consider two different 
input densities: one is flat corresponding to the naive Regge prediction and the 
other is steep in the small x region. The coefficient function can not be included 
consistently at present since the anomalous dimension has been calculated only at 
the LL order. However we consider the effects of the coefficient function because 
we could firstly clarify the above controversial aspect and secondly get some idea 
about the magnitude of the NLL order corrections in the resummation approach. 
Details of the calculations may be found in ref. [Q. 

2 Resummation of In x terms 

The flavor non-singlet part of gi in the moment space is given by, 

g,{Q\N) ^ f'dxx'^-'g.iQ^x) = ^-^C{a,{Q'), N)Aq{Q\ N) , 
Jo I 
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where 

Jo Jo 

and Aq{Q'^,x) (C) is the flavor non-singlet combination of the polarized parton 
densities (the coefficient function), (e^) is the average of quark's electric charge. 
The DGLAP equation reads, 

Q2^Ag(Q^ N) = -7(«.(Q^), N)AqiQ\ N) . (1) 

Here the anomalous dimension 7 is the moment of the "splitting" function. 

The coefficient function C{as, N) and the anomalous dimension 7(0^, A^) may 
be expanded in the powers of 



'■SI 

00 



C{a.,N) = l + Y.c\N)a>; , ^{a., N) = Y^^^N) 

k=l k=l 



a': 



where cxg = . The singular behaviors of these functions as x appear as the 
pole singularities at = 0. When x is finite, it may be enough to compute them 
to the fixed-order of perturbation. In the small x region, however, the fixed-order 
calculation becomes questionable since there appear ln"x type corrections. If 
these ln"x terms compensate the smallness of the coupling constant we must 
resum the perturbative series to the all orders to get a reliable prediction. The 
explicit next-to-leading order (NLO) calculations |]rT| in the MS scheme show 
a strong singularity at A^ = 0, 

^\N) = ^SCI-2CaCj.)^^+o(^] . (2) 



ATS 

These strong singularities (double logarithmic corrections) will persist to all or- 
ders of perturbative series. Indeed, at the k-th loop, the anomalous dimension 
and the coefficient function are expected to behave as, 

' (3) 

Our task is to resum these terms to all-orders in the perturbative expansion. 

The resummation of Inx singularities for gi has been done in Refs. ^ 0. 
The result for the "parton (quark) " target with the fixed aj^ reads, 

p2 .c+ioo fjAT / Q2\foWI^-"^ AT 



2 Jc-^oo 27ri J N-fQ{N)/8n' 



tin the genuine LL approximation, the strong couphng constant should be taken as a fixed 
parameter. 
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where /i is an arbitrary mass scale which regularizes the infrared and/or mass 
singularities. From this expression we could identify the resummed anomalous 
dimension 7 and the coefficient function C to be, 



7(«s,iV) 



hm7KiV) = -/o-(iV)/87r^ 



lim C(a„ A^) = — 



N 



Here /q is given by, 



U{N)=4n'N\l-Jl-8CF 



as 
Ar2 



1 - 



1 



ft\N) 



(4) 
(5) 



with 



N 



Dp{z) is the parabolic cylinder function. 

Now it will be instructive to re-expand Eqs. (|,^ in terms of to see whether 
these formulae sum up the most singular terms of the perturbative series. The 
expressions expanded up to 0{a^) read. 
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2Cf ]+ACf{Cf + 



C 




1 +2Cf 



k=l 



k=0 



(7) 



These results coincide with the previous expectation of Eq.(^. Furthermore one 
can easily see that the resummed expressions Eqs.(H,^ reproduce the known NLO 
results Eqs. in the MS scheme |. Therefore, it is quite plausible that Eqs. (^,|^) 
correctly sum up the "leading" singularities to all orders. 

Here a comment is in order concerning the scheme dependence. The anoma- 
lous dimension and the coefficient function individually depend on the factoriza- 
tion scheme. (Unfortunately we do not have by now any appropriate factorization 
theorems to the problem discussed in this report.) In particular, the resummed 
"coefficient function" does not have any physical meaning until the scheme de- 
pendent part of the anomalous dimension is calculated in the same scheme. To 

■fThis statement seems questionable beyond the NLO as stressed by J. Bliimlein. 
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clarify this issue, let us write the above results in the form which corresponds to 
the so-called DIS scheme. The parton densities and anomalous dimension in the 
DIS scheme are obtained by making the transformations, 

Aq Aq^^^ = CAq , 7^^^ = C-fC~^ - J 7^ In C . 

das 

Using the resummed 7 and C Eqs.(||,|^, we get the resummed part of the anoma- 
lous dimension in the DIS scheme, 

7^^^ = AT £7^ (^)\poN'td' 

k=l ^ k=2 

where the second terms come from the resummed coefficient function and d'' 
are numerical numbers independent of A^. The above equation tells us that the 
resummed coefficient function belongs to the NLL order corrections in the context 
of the resummation approach. Then, one must include the NLL order anomalous 
dimension, which has not yet been available, to see the NLL effects. 



3 Numerical Analysis 

Now we come to the numerical analyses to show how the final results are sensitive 
to the choice of the input parton densities. In conjunction with the claim in 
Ref . , we also consider the effects from the resummed coefficient function with 
the hope that the inclusion of it could shed some light on the size of the NLL 
corrections. 

Our starting point is the expression, 

N dN ;v / f^-iQ^) das nT^\ 

Jc-ioo 2m \ JoLsiQl) P j 

The anomalous dimension 7^-'^'^ which includes the resummation of ln"x terms is 
organized as follows, 

7^^5^^) _ a,7^(iV) + ah'(iV) + ^(^V, «s) - In (l + «,c^ + E{N, , 

(10) 

where 7^'^ and are respectively the exact anomalous dimension and coefficient 
function at the one and two-loop fixed order perturbation theory. K{N, as) 
{H{N,as)) is the resummed anomalous dimension Eq.(^) (Eq.(|^) with k = 1,2 
{k = 0, 1) terms being subtracted. 

It is to be noted here that the anomalous dimension at = 1 should vanish 
due to the conservation of the (non-singlet) axial vector current. The perturbation 
theory guarantees this symmetry order by order in the as expansion. However, 
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the resummation of the leading singularities in N does not respect this symmetry. 
Therefore, we need to restore this symmetry "by hand". Our choice is 



K{N,a,)^K{N,a,){l-N) . 
Of course, this is not a unique prescription |. We have used the technique in 



Ref. |T2[ to perform the Mellin inversion Eq.(g). 

We choose the starting value of the evolution to be Ql = AGeV'^. We calculate 
the evolution for two types of the input densities A and B: A is a function 
which is flat at small x ~ 0), and B is slightly steep (a ~ —0.2). The 

explicit parametrization we use is @], 

Ag(Q2, x) = N{a, p, a)r/x"(l - x f{l + ax) , 

where is a normalization factor such that / dxNx°'{l — x)^{l + ax) = 1 and r] = 
^dA {qa = 1-26) in accordance with the Bjorken sum rule. A and B correspond 
to the following values of parameters, 

A{B) : a = +0.0 (-0.2) , (3 = 3.09 (3.15) , a = 2.23 (2.72) . 

We put the flavor number Uf = 4 and Aqcd = 0.23GeV. 

First we estimate the case which includes only the LL correction 7. The evolu- 
tion kernel in this case is obtained by dropping H{N, as) in Eq. (p!0|) . Fig. la (lb) 
shows the LL results (dashed curves) after evolving to = 10, 10^, lO'^GeV^ 
from the A (B) input density (dot-dashed line). The solid curves are the pre- 
dictions of the NLO-DGLAP evolution. These results show a tiny enhancement 
compared with the NLO-DGLAP analysis and are consistent with those in Ref. 
The enhancement is, as expected, bigger when the input density is flatter. How- 
ever any significant differences are not seen between the results from different 
input densities. 

Next, we include the NLL corrections from the resummed "coefficient func- 
tion". We show the results in Fig. 2 by the dashed curves. The results are rather 
surprising. The inclusion of the coefficient function leads to a strong suppression 
on the evolution of the structure function at small x. Since the effects from the 
coefficient function fall in the NLL level, the LL terms are expected to (should) 
dominate at the small x. However our results imply that the LL approxima- 
tion is not sensible in the small x region we are interested in. As the resummed 
coefficient function is only a part of the NLL correction, we can not present a 
definite conclusion on the (full) NLL correction. But it is obvious that the NLL 
correction is very important in the experimentally accessible region of x. 

§Our final conclusion remains the same qualitatively if we choose other prescription. 
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Figure 1: The LL evolution as compared to the DGLAP results with the flat 
input A (la) and steep one B (lb). 




Figure 2: The NLL evolution as compared to the DGLAP results with the fiat 
input A (2a) and steep one B (2b). 
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4 Discussion 



In the previous section, we have shown that although the LL resummed effect 
is very small in the experimentally accessible region of x, a part of the NLL 
resummed contribution from the coefficient function drastically changes the pre- 
dictions. To understand these numerical results, it will be helpful to remember 
the perturbative expansion of the resummed results Eqs.(||,|^. Using the explicit 
values A''^; = 3,Cf = 4/3, we obtain for the anomalous dimension in the DIS 
scheme Eq.(P|), 
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DIS 



N 



-0.212 



iV2 



-0.068 ( - 0.017 (^;^ 



0.029 (^)V 



0.1411^) + 0.119 (^) + 0.069 (^) + 



(11) 



Here note that: (1) the perturbative coefficients of the LL terms (the first part of 
Eq. ([Tl|) ) are negative and those of the higher orders are rather small number. This 
implies that the LL corrections push up the structure function compared to the 
fixed-order DGLAP evolution, but the deviations are expected to be small. (2) 
the perturbative ones from the NLL terms (the second part of Eq.(|rT])), however, 
are positive and somehow large compared with those of the LL terms. This 
positivity of the NLL terms has the effect of decreasing the structure function. 
This fact that the coefficients with both sign appear in the anomalous dimension 
should be contrasted with the case of the unpolarized structure function p. 

Now it might be also helpful to assume that the saddle-point dominates the 
Mellin inversion Eq.(^. We have numerically estimated the approximate position 
of the saddle-point and found that the saddle-point stays around A^^gp ~ 0.31 in 
the region of x ~ 10^^ to 10^^. By looking at the explicit values of the coefficients 
in Eq.(|l^), the position of the saddle-point seems to suggest that the NLL terms 
can not be neglected. Since the coefficients from the higher order terms are not so 
large numerically, it is also expected that the terms which lead to sizable effects on 
the evolution may be only first few terms in the perturbative series in the region 
of X we are interested in. We have checked that the inclusion of the first few 
terms in Eq.(0) already reproduces the results of section 3. Fig. 3a (3b) shows 
the numerical results of the contribution from each terms of the NLL corrections 
in Eq.(|n]) at = lO^Ge^^ with the A (B) type input density. The solid (dot- 
dashed) line corresponds to the NLL (LL) result. The long-dashed, dashed and 
dotted lines correspond respectively to the case in which the terms up to the 
order a^, af, a^, are kept in the NLL contributions. One can see that the dotted 
line already coincides with the full NLL (solid) line. These considerations could 
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help us to understand why the NLL corrections turns out to give large effects. 
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Figure 3: Contributions from the fixed order terms in the NLL resummation with 
the flat input A (3a) and steep one B (3b). 



The final discussion concerns the convergence issue of the perturbative series 
As discussed in Refs 
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one must be careful when applying the perturbative 
approach to small x evolution. So according to Refs. we have tried to solve 
the evolution in x space with first several terms of the perturbative expansion 
being kept and what we found is that the conclusion does not change. The 
numerical results remain essentially the same. 

In summary, we have performed numerical studies for the flavor non-singlet 
Qi at small x by incorporating the all-order resummed anomalous dimension and 
a part of the NLL corrections from the resummed coefficient function. Including 
only the resummed coefficient part is not theoretically consistent, and so one 
should take into account also the anomalous dimension at the NLL level. How- 
ever, our results suggest that the LL analysis is unstable, in the sense that a large 
suppression effect comes from the resummed coefficient function which should be 
NLL correction. We need a full NLL analysis to make a definite conclusion. 
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